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Abstract 

Many observations suggest that much of the matter of the universe is non- 
baryonic. Recently, the DAMA Nal dark matter direct detection experiment re- 
ported an annual modulation in their event rate consistent with a WIMP relic. 
However, the Cryogenic Dark Matter Search (CDMS) Ge experiment excludes most 
of the region preferred by DAMA. We demonstrate that if the dark matter can 
only scatter by making a transition to a slightly heavier state (Am ~ lOOkeV), the 
experiments are no longer in conflict. Moreover, differences in the energy spectrum 
of nuclear recoil events could distinguish such a scenario from the standard WIMP 
scenario. Finally, we discuss the sneutrino as a candidate for inelastic dark matter 
in supersymmetric theories. 
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1 Introduction 



A central task of modern cosmology is to determine what the universe made of. A number 
of observations suggest that the bulk of the matter in the universe is not luminous IJ. 
Direct searches for baryonic matter in the form of massive compact halo objects (MA- 
CHOs) cannot account for the matter that seems necessary to explain these observations 

§• 

An alternative explanation is that weakly interacting massive particles (WIMPs) exist 
copiously in the halo of our galaxy but only rarely interact with ordinary matter ||. 
Candidate WIMPs from particle theory include the axion and the lightest supersymmetric 
particle (LSP) in supersymmetric theories with R-parity conservation. 

Numerous experiments have been set up in attempts to directly detect WIMPs [§], ||, §| . 
The two which are sensitive to the smallest spin independent cross sections are the CDMS 
Ge experiment |j| and the DAMA Nal experiment |J. Recently, DAMA reported the 
presence of a signal consistent with a WIMP at better than 4a. When interpreted as a 
standard WIMP with spin independent interactions, CDMS rules out nearly all of the 
DAMA 3a preferred region at 90% confidence and all of it at 84% confidence. Attempts 
to reconcile these experiments using spin dependent interactions have been shown to be 
in gross conflict with indirect detection experiments and previous direct searches 0. 

In this paper, we will show that a simple modification to the properties of the dark 
matter particle can change the kinematics of the scattering sufficiently to reconcile the 
two experiments. In particular, we explore the possibility of inelastic dark matter: relic 
particles that cannot scatter elastically off of nuclei. The outline of the paper is as follows: 
we begin by comparing the details of the two experiments and give a naive argument why 
inelastic dark matter can reconcile them. In section two we explicitly calculate the event 
rate at CDMS and DAMA taking into account the inelasticity of the scattering. In section 
three we use this calculation to study what differences can arise relative to the elastic case 
and to examine whether there are regions of parameter space that give a signal at DAMA 
but a null result at CDMS. In section four we discuss how inelastic dark matter could 
arise from a massive complex scalar split into two approximately degenerate real scalars, or 
from a Dirac fermion split into two approximately degenerate Majorana fermions. We also 
present a specific model, featuring a real component of the sneutrino as the dark matter, 
in which the mass splitting required to reconcile DAMA and CDMS arises naturally. In 
section five we discuss direct detection possibilities at future experiments. 

1.1 CDMS and DAMA 

If we are to understand the DAMA signal as evidence of dark matter, but simultaneously 
accept the null result of CDMS, we must reconsider some basic element of the WIMP 
hypothesis. Before we address such a modification, we should understand the differences 
between the DAMA and CDMS experiments. 

The DAMA experiment utilizes a set of Nal crystals at the Gran Sasso National 
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Laboratory of INFN to search for WIMPs. The basic premise of the experiment is that if 
WIMPs are present in the galaxy, as the galaxy rotates we feel a "wind" of WIMPs which 
will scatter elastically off of the target nuclei. As the Earth moves in its orbit about the 
sun, the flux and velocity distribution (as seen by a terrestrial observer) vary Rather than 
attempt to directly discriminate signal events against background, the DAMA experiment 
seeks to measure this modulation. There are two basic controls to this experiment. First, 
the signal phase must coincide with the Earth's motion in the solar system, which moves 
maximally with the galactic rotation on June 2, and maximally against on December 2. 
The second requirement is that the signal must lie dominantly in the lowest energy bins 
- a characteristic signal of WIMP scattering. 

In contrast, CDMS uses a smaller Ge target, but has excellent background rejection 
capable of distinguishing nuclear recoils from electron scatterings for scattering energies 
greater than lOkeV. As a consequence their limits are comparable to those that would 
have been expected from a null DAMA result. 

Exclusion plots are typically given in the m x — o n plane, where m x is the mass of the 
candidate and o n is the scattering cross section per nucleon. Implicit is the assumption 
that there are no great modifications in the scattering process between the two experi- 
ments. 

However, if the dark matter cannot scatter elastically, then kinematical effects sub- 
stantially distinguish the experiments. Consider two states, X- an d X+i with x+ only 
slightly heavier than x~i such that X- can only scatter by transitioning to x+- It is a 
simple kinematical constraint that X- can only scatter inelastically off of a nucleus with 
mass rriN if 

s < f3 2 ™ x m N 
2(m x + m N ) ' 

where 5 is the mass splitting between x~ an d X+- The possibility of evading direct 
detection by having a large enough splitting 5 was pointed out in [|]. Here we focus on 
the fact the constraint of equation (P becomes increasingly severe as is decreased. 
Since iodine has an atomic number of 127, while germanium has an atomic number of 
73, we have the prospect of a situation where particles will scatter at DAMA but not at 
CDMS. For /3c ~ 220km/s (a typical dark matter particle velocity), and m x = lOOGeV, 
the limits are llkeV for CDMS and 15keV for DAMA. If the mass splitting 5 were 13keV, 
such a particle would be visible to DAMA but not CDMS. 

Of course, in the halo of the galaxy there is a distribution of velocities, so the calcula- 
tion is not as simple as we have just illustrated. In the full calculation, we will find that 
the values of 5 relevant for reconciling the experiments are somewhat larger than 15 keV, 
and that the window for 5 has a size ~ 50 — 100 keV rather than ~ 5 keV. 
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2 Direct Detection Rates 



In this section we review the standard calculation of event rates at direct detection ex- 
periments |J. The differential rate per unit detector mass is given by 

dR „ T p Y f , „, . da 

Here E R is the recoil energy of the target nucleus, is the number of target nuclei 
per unit mass, p x is the local density of dark matter particles of mass m x , is the 
differential cross section for relic-nucleus scattering, and v and f(v) are the relic speed 
and speed distribution function in the detector rest frame. We take p x = .3 GeV/cm 3 . 

Because we are interested in spin-independent scattering, the differential cross section 
may be written 

da m N a n (f p Z + f n (A - Z)f 2 



dE R 2v* ni P n 

where m N is the nucleus mass, \i n is the reduced mass of the relic-nucleon system, f n and 
f p are the relative coupling strengths to neutrons and protons, and a n is the relic-neutron 
cross section at zero momentum transfer, in the elastic (5 = 0) limit. We use the Helm 
form factor [|l(| 

m)= fwy e ^ (4) 



with q = \j2m^ER, s — 1 fm, r = \/r 2 — 5s 2 , and r = 1.2A 3 . 

We assume a standard Maxwell-Boltzman distribution for the relic velocities in the 
galactic rest frame, with a root-mean-squared velocity v rms = yf^o, where we take 

Vq = 220 km/s to be the rotational speed of the local standard of rest (LSR). In our 
calculation we take the escape velocity to infinity for simplicity, when one really should 
take v esc ~ 650 km/s. By doing so we overestimate the signal for large values of the mass 
splitting 5. For a 100 GeV relic, this is a 10% effect at CDMS for 5 = 100 keV and a 
factor of two effect for S = 150 keV. Because iodine is heavier than germanium, the effect 
is far milder at DAMA, roughly 10% at 5 = 150 keV. 
The earth's speed relative to the galactic rest frame is 

V e = Vq + f or .fe COS 7 COS ((j(t — t )) . (5) 

Here v Q = v + 12 km/s, v or b = 30 km/s, u = 27r/year, t — June 2nd, and cos 7 = 
.51. Defining the dimensionless variables rj = v e /v and x min = v min /v , performing the 
velocity integration in (Q), and applying (|3[), one obtains 

dR N T m N p x 2 a n (f p Z + f n (A - Z)f ( erf(x min +77) - erf(x mi „ - 77) \ 

-b (E R ) — . (b) 



dE R Av m x p? n fl \ rj 

For the DAMA detector, one should take into account there being two species of target 
nuclei with different quenching factors. 
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Often one considers the case where f n = f p (so that the rate is proportional to A 2 ), 
and presents results in the m x - a n plane. Below we will be particularly interested in 
models in which the scattering is dominated by vector interactions arising from Z boson 
exchange, giving f n / f p = —(1 — 4 sin 2 6^) ~ —.08 (and yielding a rate that is instead 
nearly proportional to (A — Z) 2 ). In all of our calculations we take this value for /„/ f p . 

The differential rate of equation (|6]) depends on the mass splitting parameter 5 through 
Xmin, which is given by 

Xmin = — \ I — — — — + 8 J , (7) 

v Q V 2m N E R \ /i J 

where /i is the reduced mass of the relic-nucleus system. A non-zero 5 increases the 
minimum relic speed required to produce a given nuclear recoil energy. In the following 
section we explore potential consequences for direct detection signals arising due to this 
modification. 



3 Signals at CDMS and DAMA 

Before we study whether there are regions of parameter space that are consistent with both 
DAMA and CDMS, it is worthwhile to investigate the differences arising when compared 
with the elastic case. We have seen that for a given velocity of dark matter particle, it 
might be that only DAMA is be able to detect the particle, and not CDMS. Given the 
distribution of relic velocities, we can now determine what effect the inelasticity has on 
the full signal integrated over all velocities. The simplest quantity to consider is the level 
at which the signal[| is suppressed when compared with the elastic case. We plot these 
suppressions for CDMS and DAMA in figure |1|. 

We can easily see that our basic intuition is borne out. The greater the splitting 
between x~ an d X+> the greater the suppression for CDMS compared to that of DAMA. 
Since the CDMS excluded region only just covers the DAMA preferred region, even a 
factor of a few can dramatically improve the consistency of the experiments. 

However, the relative suppression is not the only relevant quantity because DAMA 
is not sensitive to the total flux, but rather to the modulation of the flux. Because 
of the inelasticity, DAMA only sees those particles on the high tail of the Maxwellian 
distribution. Consequently, a small modulation in the average velocity can lead to much 
higher modulation for a given signal when compared with the elastic case. This effect is 
demonstrated in figure 0. The combination of these two effects results in DAMA having 
significant regions of sensitivity that are inaccessible to the existing Ge experiments. 

The DAMA signal can be decomposed into background, unmodulated signal and mod- 
ulated signal as 

fi k = h + So.fc + Smj, cos(wt), (8) 

1 For our purposes here, we will consider the signal to be the events falling in the lOkeV— lOOkeV bins 
for CDMS and 2keV - lOkeV for DAMA. 
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Figure 1: Ratio of total events in iWIMP scenario to ordinary WIMP as a function of 
splitting 5 for DAMA (solid line) and CDMS (dashed), with m x = 50GeV. For DAM A 
we have integrated the total events in the 2 — lOkeV energy region, while for CDMS we 
have integrated in the 10 — lOOkeV region. For large 5 (> 100 keV), the finite value of 
the the galactic escape velocity can become important, yielding larger suppressions than 
shown. This effect is stronger for CDMS than for DAMA. 





Figure 2: Annual modulation of event rate with average normalized to one in the inelastic 
WIMP scenario (solid line) and standard WIMP scenario (dashed), with 5 = lOOkeV and 
m x = 50GeV. 
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Figure 3: Normalized modulation (S m ) as a function of energy for ordinary WIMP scenario 
(solid), inelastic WIMP scenario with 5 = lOOkeV (dashed), and inelastic WIMP scenario 
with S = 150keV (dotted), all with m x = 60GeV. 

where k indexes the energy bin of each piece of the total measured events fi. Recently, 
the DAMA collaboration published its best fit values of So,fc and S m ^ for the energy bins 
2 — 3keV, 3 — 4keV, 4 — 5keV and 5 — 6keV ||. It is tempting to fit the inelastic scattering 
case to these values, but to do so would be misleading. These best fit values are derived 
assuming the energy spectrum and relative size of the modulated piece to be given by the 
known relations for an elastically scattering WIMP. We have already seen in figure [| that 
the standard WIMP and inelastic WIMP cases can lead to very different predictions for 
the relative size of the modulated piece, so any fit to the published best-fit values would 
not be rigorous. 

Another, potentially more significant reason that we cannot use the standard WIMP 
S m ,k values comes from changes in the energy spectrum of the events. Because the scat- 
tering is inelastic, the total number of events may not rise exponentially at low energy. 
In other cases, the spectra will be nearly identical. As examples we compare in figure ^| 
the expected WIMP spectrum of the modulation signal to the spectrum in the inelastic 
WIMP scenario for two values of 5. The potential differences revealed in figure |3] make 
it possible to fit only to the model independent data recently published. We will discuss 
the details of this fit shortly. 

These spectrum differences carry over to germanium experiments. As we show in 
figure |], the changes can again be significant, and can again alter the interpretation of 
the experimental data. For instance, in the elastic case one expects an exponential rise 
in the number of relic scattering events for lower energies. Were CDMS to see many 
events in the 40 — 60keV bin, but essentially an absence of events below 40keV, this 
would be inconsistent with an elastic dark matter signal, but not with an inelastic dark 
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Figure 4: Normalized spectrum of events at CDMS for ordinary WIMP (solid) and in- 
elastic WIMP (dashed) with 5 = lOOkeV, both with m x = 50GeV. 



matter signal. Again, we do not perform a rigorous fit to the CDMS data as this would 
require an ability to accurately simulate the correlation of multiple scatterers with single 
scatterers, which we lack. 

For the purposes of generating allowed regions, we will thus use the following limits: 
for DAMA, we will use the published model independent modulation in the 2 — 6keV 
bins of 0.088 ± 0.02counts/day/kg [BJ and consider the three sigma region to be allowed. 
DAMA claims not to have modulation in the higher energy bins. Although the measured 
modulation for energies above 6keV is not published, we will take an upper limit of 
0.003counts/day/kg, which we consider quite reasonable given the errors on the best fit 
values for the higher energy bins. For CDMS, we will require a predicted mean of fewer 
than six total events, consistent with the published limits ||. 

DAMA has also reported null results arising from a pulse shape analysis (PSA) of a 
portion of their Nal data [T3| and of data from an experiment with Xe (A = 129) [[HJ]. 
Using the pulse shape, they can discriminate signal from background, and place a limit 
on the total number of events. Both of these studies affect the elastic WIMP preferred 



region for DAMA. Extracting rate limits from tables and plots of [|T^] and [jT4|, we find 
the Xe studies have the dominant impact on our allowed regions. For the Xe experiment 
we require the signal to be less than 0.7 counts/day/kg for the 13 — 15keV bin, 0.25 
counts/day/kg for the 15 — 20keV bin, 0.15 counts/day/kg for the 20 — 25keV bin, and 



0.075 counts/day/kg for the 25 — 30keV bin, consistent with published limits |H |. 

We show the allowed regions subject to these constraints for various values of m x in 
figure ||. As expected, there are broad regions that fit the DAMA data and which are 
not excluded by CDMS. It is important to note that our qualitative results are not very 
sensitive to the details of the criteria used to determine what signals are consistent with 
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Energy 


iWIMP 


DAMA 


2 - 3keV 


0.021 


0.023 ±0.006 


3 - 4keV 


0.014 


0.013 ±0.002 


4 - 5keV 


0.007 


0.007 ±0.001 


5 - 6keV 


0.003 


0.003 ±0.001 



Table 1: Binned signal rates for an inelastic WIMP with m x = 70GeV, 5 = 105keV and 
a n = 5 x 10~ 40 cm 2 , compared with the DAMA best fit values for a standard WIMP. 
CDMS would have seen an expected 0.5 events. 



experiment. The general features of figure ||] remain essentially intact even if we are more 
conservative in our estimates of the allowed counts at CDMS, or of the accuracy of the 
measured modulation at DAMA. 

As an explicit example, let us consider the point m x = 70GeV, 5 = 105keV and 
a n = 5 x 10~ 40 cm 2 . Here the modulation is quite consistent with the DAMA best fit point, 
but CDMS has only an expected signal of 0.5 events, and the Xe pulse shape analysis 
constraints are evaded. A comparison between the inelastic point and the DAMA best fit 
values for the elastic case is given in table [I]. 



3.1 Cosmological Uncertainties 

Unlike the ordinary WIMP scenario, the only inelastic WIMPs that scatter in existing 
experiments are those on the high end of the Maxwell-Boltzman velocity distribution. As 
such, there is greater uncertainty in the precise values of cr„ that fit the data than for an 
ordinary WIMP. 

In particular, there is significant uncertainty in the dispersion velocity v rms and in 
the local halo velocity v$. Although these uncertainties are O(10%), the effects can 
be amplified because of the presence of the exponential in the distribution. We have 
investigated these effects and found that the preferred cross sections can can shift by as 
much as a factor of three for m x = lOOGeV and a factor of seven for m x = 50GeV. 
Likewise the local density p x is uncertain to a factor of approximately two, and moreover, 
the presence of substructure in the halo of the galaxy can lead to amplifications of the 
local density relative to the average halo density by a factor of three or more [|ll[ . 

Finally, we must restate that we have assumed a Maxwell-Boltzmann distribution, 
which arises in the isothermal sphere model of dark matter. Changes to the velocity profile 
of the dark matter can have significant effects on the modulation for standard WIMPs 



12| , and inelastic dark matter potentially is even more sensitive to these changes. Such 
uncertainties are difficult to quanitify and we do not discuss them further. 

Altogether these uncertainties can amount to a change in the preferred o n values, but 
we should emphasize that the sizes of the regions of parameter space that yield consistency 
between DAMA and CDMS do not change dramatically. 
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Figure 5: Regions satisfying both DAMA and CDMS constraints in the 5 — a n plane, for 
(a) m x = 50GeV, (b) m x = lOOGeV, (c) m x = 300GeV. In each plot, the shaded region 
has an integrated signal in the 2 — 6keV energy range consistent with the DAMA 3cr 
region. The solid line gives the CDMS constraint and the dashed line gives the limit from 
an assumption of the absence of signal in the high energy bins at DAMA. The dot-dashed 
line gives the upper bound arising from Xe pulse shape analysis limits. The dark shaded 
region satisfies all constraints simultaneously. 
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4 Models of Inelastic Dark Matter 



Up to this point, we have considered inelastic dark matter as an interesting phenomeno- 
logical possibility, but have not addressed how such particles might arise in a reasonable 
model. One possibility is that the relic particle is a real scalar, so that its vector cou- 
pling to nuclei is forbidden by Bose symmetry. Consider a complex scalar (f> = -^(a + ib) 
coupled to an Abelian gauge field A^. Its vector interaction comes from 

|ZVI 2 3 -gA^ad^b - b&*a). (9) 

That is, the real scalars a and b couple to each other, but neither couples to itself. 

These real scalars are degenerate if the only mass term is — m 2 \(j)\ 2 , but introducing 
a small additional mass term — A 2 2 + h.c, splits this degeneracy^. If m is roughly 100 
GeV, and we want a splitting ~ 100 keV, then we require A 2 ~ (100 MeV) 2 . In the 
model of section [4.1| , which features a real component of a sneutrino as the dark matter, 
this scale for A arises naturally. 

Before discussing this model, we note that the inelastic dark matter could instead 
be fermionic. Consider a Dirac fermion ij) = ( 77 £ J that has vector and axial-vector 
couplings to quarks: 

+ g'Ai5)^qY(gv + 9Ais)q- (10) 

Assuming for simplicity that the various g's are of comparable size, the largest contribution 
to the low-energy scattering of ip off of nuclei will come from the vector- vector piece, which 
will yield an amplitude that scales roughly as the number of nucleons. The axial-axial 
piece yields a smaller spin-dependent contribution that lacks this enhancement, while the 
vector- axial pieces vanish in the extreme non-relativistic limit. 

Now suppose that in addition to a Dirac mass ~ 100 GeV for tp, the Lagrangian also 
contains a very small Majorana mass term §(?ro + f] fj), with 8 ~ 100 keV. Then the 
Majorana fermion mass eigenstates are 

mi = m — 5 (11) 

vri2 = m + 5. (12) 

The vector current essentially couples xi to X2, with only a small additional piece ~ 5/m 
coupling each mass eigenstate to itself: 

V>7mV> - i(Xi0»X2 ~ X2°vXi) + ^—(X2 W i^X2 ~ Xi^ixXi)- (13) 

Because S/m ~ 10~ 6 , we ignore the second term, and find that the only way for xi t° 
scatter coherently off of nuclei is to make a transition into the heavier X2 state. This 

2 Of course, A violates gauge invariance, and can only arise once the gauge symmetry of the theory 
has been broken. 



xi - 

X2 ^ -fiiv + O 
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inelastic process will dominate relative to the elastic, spin dependent scattering provided 
that the coherence enhancement, which gives a factor ~ A 2 ~ 5 ■ 10 3 in the cross section, 
overcomes than the suppression due to the inelasticity. In this case, the rate can depend 
sensitively on the mass of the target nucleus, as desired. 



4.1 Sneutrino Dark Matter 

Interestingly enough, a suitable candidate for inelastic dark matter has already been 
discussed in the literature. In supersymmetric theories with lepton number violation, the 
LSP can be a real component of the sneutrino []8], |I5, IBfl: a Lagrangian term — A 2 z>z> + /i.c. 
lifts the degeneracy between the sneutrino's odd and even CP eigenstates z>_ and z> + . 
This splitting prevents elastic scattering of the lightest state, z>_, off of nuclei through 
Z exchange^. There is still the challenge of achieving a cosmologically interesting relic 
abundance - for an ordinary 100 GeV sneutrino Qt comes out too small. In ||, this 
problem was resolved by taking the mass splitting between z>_ and + to be large enough 
to prevent coannihilation via an s-channel Z in the early universe, 5 > 5 GeV, leading to 
a radiatively generated neutrino neutrino mass m v > 5 MeV. Different approaches were 



taken in the models of |L5], 16, 18|. These models feature standard model singlet scalars 



n that are kept light by a global symmetry [15, 16] in analogue to the Giudice-Masiero 



solution to the \i problem |L7j, or by a gauged B-L symmetry [fLql . The singlet states 
mix with ordinary sneutrinos through weak scale A terms, so that the gauge interactions 
of the mass eigenstates are suppressed by mixing angles. This suppression allows for 
an interesting relic abundance even for values of 5 too small to prevent coannihilation 
between z>_ and u, . 



For concreteness we will specialize to the model of | lq| . The global symmetry that 
prevents a tree level mass for the singlet n states is broken by the vev of a spurion X 
that also breaks supersymmetry. We assume that the A and F components of X both 
have intermediate scale vevs: (Ax) ~ a/ (Fx) ~ rrij ~ y/vMpi. The spurion couples to 
the neutrino and singlet superfields according to 

i i / v-f v \ 

C D — [XLNH U ] F + — [XtjViV + + + h .c. (14) 

The operators of (|14| ) can be justified by ordinary R parity (under which N is odd and 
X is even), together with an R symmetry where N has R charge 2/3, X has charge 4/3, 
and L and H u have R charge 0. As discussed in |TEfl , at tree level QT4] ) yields a neutrino 
mass matrix whose light eigenvalue is ~ v 2 jMp\. However, ([14]) also contains 

C D -Alnh u - A 2 (nn + h.c), (15) 



3 There are contributions that will induce an elastic scattering, for instance from Higgs exchange, but 
these are all small and can be ignored for our purposes here. 
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with A roughly weak scale and A 2 ~ mf/Mp>j. These interactions radiatively induce a 
Majorana mass for the left-handed neutrino 

g 2 v 3 / 2 , . 

m u ~ Q0 . 2 1/2 (16) 
384vr 2 Mpf 

that is larger than that obtained from the tree-level seesaw, and moreover, in roughly the 
correct range for explaining the atmospheric neutrino anomaly. 

For our present purpose, however, the impact of flTSP on the scalar masses is what 
matters most. Neglecting the small lepton number violating mass parameter A, the 
sneutrino mass-squared matrix is 

/ \ ( m2 T -ha Av sin 8 \ ( i> \ , s 

v ' \ -^Avsm8 m l R I \ n * J v 1 

The A term coupling induces a mixing between v and n, yielding a lighter mass eigenstate 

v\ = — z>sin# + n* cos9. (18) 

The coupling of U\ to the Z boson is thus suppressed by sin 2 9. The lepton number 
violating parameter A 2 lifts the degeneracy between the CP-even and CP-odd components 
of Ui, leading to a small mass difference 

A 2 

5~2cos 2 # — . (19) 
mi 

For this splitting to resolve the conflict between CDMS and DAMA, one needs 5 ~ 50—100 
keV, roughlyQ. For a 100 GeV sneutrino, this implies [X^XX% ~ mf ~ (3- 10 10 GeV) 5 , 
corresponding to a reasonable value for the intermediate scale. 

To explore the feasibility of this scenario, we apply the same criteria used in section 
3 to establish consistency with CDMS and DAMA for 5 = 50 and 100 keV, and display 
the allowed regions in the (m?, sin9) plane. Note that because the scattering off of nuclei 
is suppressed both by the inelasticity of the reaction and by a sin 4 9 factor, the ability 
to obtain a large enough signal at DAMA depends crucially on the fact that ordinary 
sneutrinos give a signal roughly three orders of magnitude above present bounds. We also 
calculate the relic abundance as a function of and sin 9 using standard methods. The 
results shown in figures |§a)-d) indicate that there are indeed regions of parameter space 
featuring interesting relic abundances and acceptable direct detection rates. In the early 
universe, the efficiency of annihilation processes that occur via s-channel Higgs exchange, 
such as z>_P_ — ► bb, ZZ, W + W~ , is sensitive to the size of the trilinear scalar coupling 
A, leading to the dependence of the relic abundance on A evident in the figures. 

4 The lifetime of v + is r ~ (jp^") ( 100 ^ kcV ) J (4 . 1Q 2 years), so for the mass splittings and mixing 



s 

angles of interest, it is safe to assume that only P_ is present today. Photons can be produced in these 
decays, but the decays take place before recombination for the parameters of interest, and the photons 
are soft enough to render negligible the effect on the CMBR spectrum. 
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a) b) 



0.6 




m~ /GeV m~ /GeV 

Figure 6: For the sneutrino dark matter case, regions that satisfy the direct detection 
requirements of section 3, plotted along with filled contours of Q^h 2 . The lighter shaded 
region corresponds to .05 < Q^h 2 < .3 and the darker shaded region corresponds to 
.3 < Q„h 2 < .64. The region between the dotted contours has an integrated signal in the 
2-6 keV range consistent with the DAMA 3cr region. The solid line gives the CDMS 
constraint, the dashed line gives the limit from the absence of signal in the high energy 
bins at DAMA, and the dot-dashed line gives the constraint arising from Xe pulse shape 
analysis data (regions below these lines are allowed). We take 5 = 50 keV for a) and b) 
and 5 = 100 keV for c) and d). For a) and c), we use A = 25 GeV while for b) and d), 
we take A = 50 GeV. For each plot we take tan/3 = 50, = 115 GeV, and a bino mass 
of 300 GeV, with the assumption of GUT unification of gaugino masses. 
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Energy /keV 


S m ,fc(cpd/kg/keV) 


DAMA best 


inelastic v 


2-3 


.027 


.027 


3-4 


.013 


.013 


4-5 


.005 


.006 


5-6 


.002 


.002 



Table 2: S m ^ values obtained using DAMA's best fit point, o = 7.2- 10 -6 pb and m x = 52 
GeV for the standard WIMP case, and values obtained taking 5 = 70 keV, = 70 GeV, 
and sin 4 # = 1/70 for the sneutrino inelastic dark matter case. 



Just as one specific illustrative example, consider the parameters mo = 70 GeV, 5 = 70 
keV, and sin 4 # = 1/70. This choice of sin# leads to an interesting relic abundance for 
a broad range of SUSY parameters. For this choice of mp, S, and sin^, we calculate a 
mean of less than 2 events at CDMS and satisfy the constraints from the Xe pulse shape 
analysis. Moreover, as shown in table 0, the values of S m ^ we obtain for DAMA in the 2-6 
keV energy range are nearly identical to those we obtain in the elastic case using DAMA's 
best fit point a n = 7.2 ■ 10~ 6 pb and m x = 52 GeV. 



4.2 Indirect detection 

As dark matter passes through the sun, it can scatter off of nuclei and be captured in the 
sun's potential well []B| . After a significant amount of dark matter has been captured, it 
can annihilate into other particles. If muon neutrinos are produced, those that reach the 
Earth can produce high energy muons through charged-current interactions. A number 
of experiments have attempted to direct WIMP matter indirectly by looking for these 



upward-going muons, leading to a current limit on their flux of 10 cm s [£0|, ^1], P2|| . 



Even within the model of section fO, the expected flux of upward-going muons is 



quite uncertain, for a number of reasons. First, the capture rate in the sun is sensitive 
to the parameters m^, sin 6 and 5. Second, if the captured sneutrinos annihilate directly 
into neutrinos, the flavor of the lightest sneutrino determines what flavor of neutrino is 
produced, and details of the neutrino masses and mixings impact the flavor of the neu- 
trino detected at the Earth. Third, cosmological uncertainties mentioned in section [TT] 
can change the preferred region of a n and thus the capture rate. Finally, relatively minor 



extensions to the model of section |4.1| can also complicate matters. As a consequence of 
these various sources of uncertainty, indirect techniques do not rule out inelastic sneutrino 
dark matter. However, they do impose strong constraints, as broad regions of parameter 
space lead to signals above experimental bounds. Moreover, indirect detection experi- 
ments offer the strong possibility of detection if the bound on the muon flux improves 
considerably p3 . 



One might expect that the same inelasticity that suppresses the signal at CDMS 
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should be even more effective in suppressing the capture rate by the sun, which is mainly 
composed of relatively light nuclei. In fact, this is typically not the case. Because particles 
passing through the sun are unusually energetic (the escape velocity at the surface of the 
sun is much larger than the average velocity of a halo particle), the inelasticity is less 
relevant in the sun than at direct detection experiments. 

In what follows, we have followed |24| in calculating solar capture rates and the induced 
muon flux, but have modified the approach to approximate the suppression of the capture 
rate due to the inelasticity (see Appendix). This suppression depends on the mass of the 
nucleus. For example, for 5 = 100 keV and = 100 GeV, we find a factor ~ 20 
suppression for scattering off of oxygen in the sun, and a factor ~ 2 suppression for 
scattering off of iron. 

We will separately consider first, the case in which the sneutrinos cannot annihilate 
into W's, and second, the case in which they can. For vrty < mw, sneutrinos in the 
sun typically annihilate dominantly to neutrinos via t-channel neutralino exchange. If 
we neglect cosmological uncertainties, we find that for values of m^, 5 and sin# that 
lead to interesting relic abundances and consistency with CDMS and DAMA, the flux of 
neutrinos produced is quite large. If these are all muon flavor, we would expect a flux 
of upward-going muons of at least ~ 6 x 10 _13 cm _1 s _1 , in conflict with experimental 
results. On the other hand, if the sneutrinos annihilate into electron neutrinos that do 
not oscillate into muon neutrinos, bounds from direct detection are evaded entirely. If 
this is the case, only direct detection experiments will be able to yield a positive signal. 

Finally, we note that there are specific parameter choices for which the dominant 
annihilation of sneutrinos in the sun is through s-channel Higgs to bb. We find that this 
allows the flux of upward-going muons to be as small as a factor of ~ 2 above current 
limits for parameters that yield an acceptable abundance and acceptable direct detection 
signals. 

With cosmological uncertainties included, more scenarios are allowed. Relatively small 
variations in v rms can accomodate factors of two, such as if the dominant annihilation is 
into bb. Direct annihilation into muon neutrinos would require a more specious conspiracy 
of errors. For instance, if the solar system were presently in an anomalously high density 
region of the galaxy arising from substructure, and if v rms were 3a above the value we 
have used, experiments could accomodate as much as one-third of the neutrinos produced 
being muon flavored. This seems quite unlikely, but is, at least in principle, still allowed. 

However, for these lighter sneutrinos, if indirect detection experiments improve by an 
order of magnitude, they will be able to probe almost all of the parameter space, even 
accounting for a broad class of cosmological uncertainties, and situations where there is 
annihilation to bb. 

For heavier sneutrinos (m^ > m\y), the dominant annihilation processes in the sun 
can easily be s-channel Higgs exchange to Ws and Z's. In this case we find that it is 
possible to reduce the expected signal at direct-detection experiments to a factor of ~ 3 
above current limits for parameters consistent with DAMA, CDMS, and Q„h 2 ~ .1. These 
heavier sneutrinos are less affected by cosmological uncertainties, but these uncertainties 
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still make it impossible to rule out this scenario. Future improvements in indirect detection 
could rule out this region of parameter space, especially as the experimental signal is less 
sensitive to the flavor of the lightest sneutrino than for the case of < mw 

Of course, this discussion applies only to the model of section 3~T , and one can consider 
modifications to the model that suppress the indirect detection signal. The premise of 
the model is that light standard model singlets are natural. Given this, if we add to that 
model another standard model singlet 77, with the same R-charge as the right handed 
neutrino N, but opposite R-parity, we expect a superpotential interaction r/NN. Then 
through t-channel r\ exchange, z/'s can annihilate to right-handed neutrinos. If these decay 
dominantly into muons or electrons (rather than tau's) and off-shell W's, we find that it 
is possible to bring the flux of upward going muons induced by the W decay products 
down to current limitsQ Future indirect experiments would still likely be able to see the 
decay products of these right handed neutrinos. 

This is just one example of a modification to the model which diminishes the signal, 
and there may be others, but such uncertainty is difficult to quantify. While indirect 
experiments offer a good opportunity to test specific models and regions of parameter 
space, there is an excellent likelihood that upcoming direct detection experiments will be 
able to determine whether inelastic dark matter is the resolution of the conflict between 
DAMA and CDMS. 



5 Future Experiments 

In the inelastic dark matter scenario, the boundaries of the DAMA preferred region are 
not far from the current limits from CDMS. Planned experiments should be able to cover 
the existing DAMA region. Most important are planned improvements to germanium 
experiments, and the CRESST experiment, which will use the heavy element tungsten. 

CDMS will soon be moving to the Soudan mine, and should be able to improve its 
limits by at least two orders of magnitude | 25| . The GENIUS Ge experiment should 



go well below that, likely allowing both to test much of the preferred regions discussed in 
section |3|. 

There is a caveat in this statement: in generating the plots of figure |5], we neglected 
to include the effect of a finite galactic escape velocity. This was a harmless simplification 
for our purposes there, because the effects at DAMA due to the finite galactic escape 
velocity are relatively minor. The effects can be much larger at CDMS. 

Recall that the requirement for scattering is 

5 < g mNm - . (20) 
2 mjv + rn x 

This constraint is particularly stringent for light candidates. For instance, with m x = 
50GeV and v esc = 650km/s equation (|20|) tells us that only for 5 < 122keV can one hope 

5 In such a scenario, the relic abundance is modified, but it is still possible to have fl„h 2 ~ 0.1. 
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to obtain any signal at all at a germanium detector (recall that the highest velocity of 
particles incident on the earth is v esc + v Q ). Thus, the higher 5 regions may not be testable 
at CDMS. 

For heavier candidates, the finite galactic escape velocity is not especially important, 
even at CDMS. With a galactic escape velocity of 650km/s, and m x = lOOGeV, the cutoff 
for 5 is 172keV. On the other hand, the galactic escape velocity is not a particularly well 
known quantity, and if instead we take v esc = 450km/s, the cutoff for 5 is only 102keV. 

These uncertainties make the CRESST experiment |27j, using tungsten, especially 
significant. Because tungsten (A = 183) is heavier than iodine (A = 127), given ade- 
quate exposure time, CRESST should cover the DAMA preferred region, irrespective of 
cosmological uncertainties. 

A very real possibility is that both germanium and tungsten experiments will have 
signals, which, when interepreted as elastic scatterings, would be inconsistent with one 
another. The most striking possibility of all is a spectrum deformation at the germanium 
detectors, as discussed in section |3|. If CDMS were to see an excess of events in the 
30 — 70keV region, but no excess below 30keV, it would be a compelling signature of this 
scenario. 

6 Conclusions 

If in fact the majority of the matter of the universe is non-baryonic, the attempt to 
determine its nature is one of the most exciting endeavors of modern cosmology. Existing 
dark matter searches have already begun to probe interesting regions of parameter space 
for candidate particles such as neutralinos and axions. 

The positive result from the DAMA experiment is difficult to a understand in terms 
of these candidates, as it is in seeming conflict with constraints arising from the CDMS 
experiment. We have seen that this conflict vanishes if we allow for the possibility that 
the dark matter particle can only scatter inelastically. 

We have shown that the sneutrino, when mixed with a singlet scalar with weak lepton 
number violation, is a viable candidate for inelastic dark matter. The regions of parameter 
space which give an interesting relic sneutrino abundance overlap with the regions which 
give a positive DAMA signal. Indirect detection experiments tightly constrain models of 
sneutrino dark matter, but do not rule them out. 

Even absent a particular model, we find it interesting that such a simple modification 
of the dark matter's properties can give remarkably different predictions, including the 
suppression of a signal at CDMS. We consider these results sufficiently interesting as to 
warrant an analysis of the full DAMA data set should the raw data become available. 
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A Appendix 

Here we describe how we approximate the suppression of the rate of relic capture by the 
sun due to the inelasticity of the scattering. Ignoring nuclear form factors, the scattering 
probability for a given relative velocity w is equally distributed between the minimum 
and maximum nuclear recoil energies AE min and AE max (these parameters depend on 
the nucleus mass, the relic mass m, w, and 8). Ordinarily, the low-energy scattering 
cross section is independent of w, but in the inelastic case there is an additional phase 
space factor a/1 — 25/ (pw 2 ), where \x is the reduced mass. Capture only occurs when 
AE > AE capture = 1/2 [mw 2 — (m + 5)v 2 sc {r)) — 8 holds. Here v esc {r) is the (position- 



dependent) escape velocity, which we approximate as [23 



Ves C {r) = vl-^{vl-vH (21) 

where v c = 1354 km/s, v s = 795 km/s, and M(r) is the mass contained within the radius 
r. The capture rate off of a given species of nuclei is then proportional to 

/ dr r 2 p(r) / dw w 3 e y/i - 25/(pw 2 ) ^ ax ^7 > ( 22 ) 

where p(r) is the mass density of the species and vq is the rotational speed of the local 
standard of rest. We calculate this factor (which does not account for form factor sup- 
pressions) in the elastic (8 = 0) and inelastic cases to estimate the suppression coming 
from the inelasticity. We then obtain capture rates by multiplying this suppression with 
the rate obtained for the elastic case using the formulae of [^4j (which do include form 
factor suppressions). 
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